Electronic state around vortex in a two-band superconductor 
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Based on the quasiclassical theory, we investigate the vortex state in a two-band superconductor 
with a small gap on a three dimensional Fermi surface and a large gap on a quasi-two dimensional 
one, as in MgB2. The field dependence of zero-energy density of states is compared for fields 
parallel and perpendicular to the ab plane, and the anisotropy of the vortex core shape is discussed 
for a parallel field. The Fermi surface geometry of two-bands, combining the effect of the normal- 
like electronic state on the small gap band at high fields, produces characteristic behavior in the 
anisotropy of c- and afo-directions. 
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I. INTRODUCTION 

Since there are many materials with a multi-band 
structure, it is important to study the contribution of 
the multi-band structure on superconductivity in or- 
der to understand the superconducting state in these 
materials^ Superconductivity of MgB2 has been well 
studied both experimentally and theoretically after its 
discovery^ Now it is recognized that MgB2 is a typical 
example of a multi-band superconductor, which has two 
bands and two superconducting gaps, i.e., a large super- 
conducting gap in the a band and a small gap in the 7r 
band, 3 . 4 . 5 . 6 . 7 

On the other hand, the electronic states around the 
vortex when applying a magnetic field can be a probe 
for the anisotropy of the superconducting gap and the 
Fermi surface structure. The field dependence of the 
density of states (DOS) N(E = 0) depends on the gap 
anisotropy. 8,9 From the field orientation dependence of 
N(E = 0) we can obtain information of the node posi- 
tion in anisotropic superconductors i& These features are 
examined experimentally by the specific heat or thermal 
conductivity measurements 

When this analysis is applied to MgB2, we detect the 
feature of two-band and Fermi surface anisotropy. The 
rapid increase in the field dependence of the DOS at low 
fields as N(E = 0) ~ fjo.23 anc j ^. ne i ar g e vortex core 
size observed by scanning tunneling microscopy (STM) 
in MgB2 are understood as characteristics of two-band 
superconductivityiiiii 4 * 1 ^ These come from the contribu- 
tion of the small gap band, whose electronic state be- 
comes normal-like state by applying a weak magnetic 
field, as explained by theoretical calculations for H \ 
Ci ^6 i i7 i i8 Therefore, characters of a two-band supercon- 
ductor are understood for H \\ c. 

As the next step, two-band characters are examined 
for H || ab, or by field-orientation dependence. In the 



field-orientation dependence, the contributions from the 
isotropic three dimensional (3D) Fermi surface as shown 
in Fig. HJa) are expected to be isotropic. On the 
other hand, the quasi-two dimensional (Q2D) Fermi sur- 
face as shown in Fig. Qlb) gives an anisotropic con- 
tribution, reflecting the coherence length ratio £ a b/£c, 
which is related to the Fermi velocity anisotropy ratio 
(bafc| 2 ) 1 ^ 2 /(|w c | 2 ) 1 ^ 2 - As the contributions of the cr-band 
with a Q2D Fermi surface and a 7r-band with a 3D Fermi 
surface are coupled in MgB 2 and the dominant contribu- 
tion changes depending on the field range, the study for 
H || ab or for field-orientation dependence is important 
in order to further extract the information of a two-band 
superconductor and the Fermi surface anisotropy. 
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FIG. 1: Schematic view of the isotropic 3D Fermi surface 
(a) and Q2D Fermi surface (b). We consider the cases of two 
magnetic field orientations, H \\ c and H || ab. 

Focusing on this point, some experimental studies 
and analyses were performed. In the specific heat and 
the thermal conductivity measurement jisLSS the field- 
dependences are compared for H \\ ab and H || c, which 
shows that only a slight difference between N(E = 0, H \\ 
c) and N(E = 0, H \\ ab) is observed at low fields, while 
at high fields there is a pronounced anisotropy. This 
result was explained due to the fact that the dominant 
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isotropic contribution to N(E = 0) at low fields is coming 
from the 7r band with a 3D Fermi surface, and high field 
anisotropic behaviors are due to the a band with a Q2D 
Fermi surface. The different temperature dependence of 
anisotropy ratios ^ H = H C 2,ab/H c2jC and 7a = A c /A afc is 
also a consequence of this multi-gap (multi-band) nature 
of superconductivity in MgB2 ^» There were some works 
studying the contribution of 3D and Q2D Fermi surface 
sheets on the vortex state by quasiclassical theory 18 i 22 i 23 
There, ff c2 -behavior is compared for H || c and H || ab, 
and the electronic states of the vortex state were dis- 
cussed for H || c. In this study, we discuss the electronic 
state of the vortex states mainly for H \\ ab. 

For H || ab, it is important to examine the vortex 
core shape, since the core shape is determined by the 
anisotropy ratio of the coherence length. For this pur- 
pose, we selfconsistently determine the profile of the 
order parameter in the vortex state. In this field- 
orientation H || ab, we expect a highly anisotropic core 
shape from the Q2D Fermi surface contribution of the a 
band. However, the vortex core image for H \\ ab ob- 
served by STM shows that the vortex core is a circular 
shape at a low fields This indicates that the anisotropic 
contribution of a Q2D Fermi surface is not clear and that 
the contribution of the 3D Fermi surface is dominant in 
the vortex core image of STM at this field range. 

The purpose of this paper is to investigate the elec- 
tronic state of the vortex lattice state mainly for H \\ ab 
using quasiclassical theory in the clean limit of a two- 
band superconductor with 3D and Q2D Fermi surfaces. 
The quasiclassical theory can be applied in all range of 
temperatures and magnetic fields. After describing our 
formulation of the quasiclassical theory in Sec. ^ we an- 
alyze the field dependence of N(E = 0) comparing two 
field orientations H \\ ab and H || c in Sec. IIIII and clarify 
the contributions of multi-gap superconductivity and the 
Fermi surface geometries. In Sec. IIVI we study the vor- 
tex core structure and the local density of states (LDOS), 
which are observed by STM, and discuss the anisotropy 
of the vortex core shape for H \\ ab. In Sec. El we show 
the quasiparticle excitation spectrum outside vortex core 
to see the behavior of the gap edge at finite fields. The 
last section is devoted to summary and discussions. 



II. FORMULATION 

We consider a simplified model of a two-band sys- 
tem with a large superconducting gap band and a 
small gap band (denoted as L-band and S-band, re- 
spectively). For simplicity, the superconducting gap at 
each band is assumed isotropic. We use the following 
model of the 3D and Q2D Fermi surfacesi 18 ' 22 ! 23 As 
shown in Fig. QJa), the S-band corresponding to the 
7r-band is assumed to have a spherical Fermi surface, 
given by E s (k F ) = (h 2 /2m)(k 2 x + k 2 y + k 2 z ) = E F 

with Fermi energy Ep = h 2 kp/2m. Fermi velocity is 
given by v(k FS ) = dE s /dk = v F o{k x ,k y ,k z )/k F = 



v F o(sin# cos (/>, sin^ sin0, cos 9) with < 9 < tt and 
< (j> < 2tt. As shown in Fig. ^b), the L-band is as- 
sumed to have a Q2D Fermi surface of cylinder-like shape 
with small ripples, given by £x(k F ) = (H 2 /2m)(k 2 + 
ky) — tcosk z — Ep and Fermi velocity v(kpi) = 
vfol (cos (j), sin (j), v z sinfc z ), where v z is small and related 
to the anisotropy ratio as l/v z ~ 'Jh ■ In our calculation, 
we set l/v z = 6 and vflo = vfo- 

Vortex structure and electronic states are calcu- 
lated by quasiclassical Eilenberger theory in the clean 
limit ,£L25i2£i21 First, the quasiclassical Green's functions 
giiuj n , k Fi , r), f(iu> n , bFj,r) and fHiu> n , k Fj , r) are calcu- 
lated in a unit cell of the vortex lattice by solving Eilen- 
berger equation 

{cj n + v(k Fj ) • [V + iA(r)]} / = A,-(r)«7, (1) 
{u n - v(k Fj ) • [V - iA(r)]} /t = A*(r) 5 (2) 

in the so-called explosion method, where g = (1 — fp) 1 ' 2 , 
Reg > 0, v(k F:( ) = v(k F:) )/u F o, j = L,S, and Matsubara 
frequency oj n — (2n + 1)ttT . In the symmetric gauge, 
A(r) = |B x r + a(r), where B = (0, 0, B) is a uni- 
form field and a(r) is related to the internal field b(r) as 
b(r) = V x a(r). The unit vectors of the vortex lattice 
are given by Ui = (a x ,0, 0), and u 2 = (^a x ,a y ,0). The 
pair potential and vector potential are selfconsistently 
calculated by the relations 

A;(r)=2T^ J2 YjyfihJnMs',*), (3) 

w„>0j',k Fj , 

J(r) = -2T£- 2 E v(k Fj -)Im ff (iw„,k Fj -,r), (4) 

where J(r) = V x V x A(r), and k = (7C(3)/18) 1 / 2 k B cs 
with Rieman's zeta function C(3). We set the energy cut- 
off lj c — 20T c o, kbcs = 30, and the DOS in normal state 
at each Fermi surface N l — Aos — 0.5AT . The integral 
J2j k F takes account of iVJy and the k F -dependent DOS 
(oc |v(k F:) )| _1 ) on the Fermi surface. The LDOS is given 
by 

N taM (E,r)= £ Nj(E,r) 

3=L,S 

= E Re 5( iw « -> E + irj, k F j,r), (5) 

where g is obtained by solving Eqs. JQ| and (01 under 
iu> n — > E + bj. We typically use 77 = 0.01. The spatial 
average of the LDOS gives DOS; 

N totai (E)= £ NjiE) = (N(E,r)) r . (6) 

j=L,S 

We assume that superconductivity in the S-band oc- 
curs by Cooper pair transfer V^s from the L-band. 
Therefore, we set the pairing interaction V$s = in the 
S-band, and use Vsl = 0.32Vll so that A^/Ag ~ 3 at 



a zero field. Throughout this paper, temperatures, ener- 
gies, lengths, magnetic fields, and DOS are, respectively, 
measured in units of T c o, TrksTco, Rq = hvpo/2itkBT c o, 
Bo = 0o/27T-Rq and No, where T c0 is the superconduct- 
ing transition temperature in the case Vls — Vss = 0. 
Our calculation is performed at T = 0.1T c o. The shape 
of the vortex lattice for H \\ c is the triangular lattice 
with a x /(2a y /y/3) = 1. For H \\ ab, the anisotropic ratio 
a x /(2a y /y/3) of the vortex lattice varies from 1.2 at low 
fields to "f H ~ 6 at high fields in MgB 2 i 24 i 28 We study two 
cases a x / (2a y / ya) = 1.5 and 6 for this field orientation. 

We selfconsistently calculate the order parameter 
Ai(r), Ag(r) and the vector potential a(r), under a given 
vortex lattice ratio a x /(2a y /y/3). We alternatively solve 
the Eilenberger equations [Eqs. Q and @] and the self- 
consistent condition [Eqs. J3J and Q], until a sufficiently 
selfconsistent solution is obtained. As an initial state 
of the calculation, we use the Abrikosov solution of the 
lowest Landau level for Ai(r), and set As(r) = and 
a(r) = 0. By this selfconsistent calculation, we can prop- 
erly estimate the size and shape of the vortex core. In 
the Abrikosov solution, vortex core size is always propor- 
tional to the inter-vortex distance. In our selfconsistent 
method, the vortex core size remains to be in the order of 
the coherence length even when the inter-vortex distance 
becomes large at low fields. This proper treatment of the 
vortex core size is necessary to correctly estimate the low 
energy excitations of the quasiparticles especially at low 
fields. The selfconsistently obtained a(r) does not give a 
significant contribution on the electronic state, since we 
consider the case of high-K. 



III. FIELD DEPENDENCE OF ZERO-ENERGY 
DENSITY OF STATES 

At low fields, low energy quasiparticles are localized 
around vortex cores in full-gap superconductors. With 
increasing field, the distribution of low energy quasi- 
particles around a vortex core overlaps with that com- 
ing from neighbor vortex cores, since inter-vortex dis- 
tance decreases^ Therefore, low energy quasiparticles 
extend outside of the vortex core. The zero-energy to- 
tal DOS N to tai(E — 0) is the spatial average of these 
low-energy quasiparticle distributions. Field dependence 
of total DOS Ntota\(E = 0) and each band contribution 
N L (E = 0) and N S (E = 0) are shown in Fig. ^a) both 
for H || ab and H \\ c. The DOS on the S-band increases 
rapidly at low fields, and saturates to normal state value 
0.5No at higher fields. The increase of total DOS at high 
fields is due to the L-band contribution. These are con- 
sistent with previous calculations»iS^2^ 

In this study, we compare the curves for H \\ ab and 
H || c. The rapidly increasing S band contribution N$ 
is almost the same for both orientations at low fields, 
reflecting the 3D isotropic Fermi surface shape. Slowly 
increasing the L-band contribution N L shows a large dif- 
ference, reflecting the anisotropy of H C 2 dominantly com- 



1 

II 

0.5 



/H//c ' 

• 

. F LULdi 


1 

^total 


1 < 
^H//ab '. 








* 


S 




TZ 






F P 

r^-® 

8r e ~^ i i 




i 




0.5 

ii 



(c) 



Mob 


totaUj^^ 















B 



1 



FIG. 2: (a) Field dependence of zero-energy DOS iVtotai = 
N L + N s (•) for H || ab (solid lines) and H \\ c (dashed 
lines). Large-gap band contributions Nl (°) and small-gap 
band contributions Ns (o) are also presented. Vortex lattice 
ratio a x /(2a y /V3) = 6 for H \\ ab and 1 for H \\ c. (b) Field 
dependence of maximum pair potential amplitudes |Ai| (o) 
and | As | (o) in the vortex lattice state for H || ab (solid lines) 
and H || c (dashed lines), (c) Field dependence of zero-energy 
DOS for a x /(2a y /V3) = 1.5 (thin lines with x) and 6 (thick 
lines) for H || ab at low fields. 



ing from the Q2D Fermi surface. As a result, A^otai is 
almost isotropic at low fields, and largely anisotropic at 
higher fields, which is consistent with experimental ob- 
servations iliLSSi 

Figure EJb) indicates that, while electronic states in 
the S-band are normal-like state at high fields, the pair 
potential As in the S-band persists up to H C 2 with a gap 
ratio |Ai|/|As| relatively unchanged, as well as that As 
persists up to T c at H = 04*^& This is easy to understand 
since superconductivity in the S-band is induced by the 
Cooper pair transfer from the L-band. As long as we 
have the pair potential |A^| in the L-band, the induced 
pair potential in S-band is non-vanishing. This behavior 
is similar both for H || c and H \\ ab. 

It is necessary to see how DOS depends on the 
anisotropic ratio a x / (2a y / for a unit cell of the vor- 
tex lattice, which is given in our calculation. We plot the 
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DOS for a x /{2a y /V3)=6 and 1.5 in Fig. 13c). At low 
fields, DOS does not significantly depend on the shape 
of the unit cell, because low-energy electrons are local- 
ized around vortex cores. Therefore, our numerical re- 
sults do not significantly depend on the delicate tuning 
of the vortex lattice ratio a x /(2a y /y3) at low fields. At 
higher fields, where the low-energy quasiparticles around 
vortices overlap each other, there appears a deviation 
depending on the vortex lattice shape. At these fields, 
however, the equilibrium vortex lattice in MgB 2 is a dis- 
torted hexagonal with a x /(2a y /^/S) ~ jh ~ 6. 28 There- 
fore, it is reasonable to study the vortex structure with 
a x /(2a y /y3) — 6 at high fields. 




FIG. 3: (a) Field dependence of DOS iVtotal (•), N L (o) 
and Ns (o) for H || ab (solid lines) and H || c (dashed lines) 
when both bands have Q2D Fermi surfaces with l/v z = 6. 
a x I (2a y I V3) = 6. (b) Field dependence of (a) is replotted as 
a function of 6B for H || c {B for H || ab), as ~ 6. 

To show an example that the field-dependence of A^otai 
is changed by the Fermi surface geometry in a two-band 
superconductor, we also calculate the case when both L- 
and S-bands have Q2D cylindrical Fermi surfaces with 
l/v z — 6, for comparison. In this case, as shown in Fig. 
|3fa), field-dependences of A^otai are anisotropic between 
H || ab and H || c both at low and high field, because 
S-band contributions at low fields are also anisotropic 
due to the Q2D Fermi surface contribution. When all 
bands have the same Fermi surface geometry as in this 
case, we can expect that ^-dependences of N to tai(E = 0) 
are roughly scaled by H C 2 anisotropy. In Fig. OJb), the 
field dependences of zero-energy DOS are replotted as a 
function of 6B for H || c (B for H || ab) as ~ 6. There, 
calculated data of AT tota j (and also Nl, Ns) for H \\ ab 
and H \\ c are on the same curve, i.e., field dependences 
are well scaled by H C 2 anisotropy. On the other hand, 
the field dependences in Fig. |2Ja) are not scaled by the 
H C 2 ratio, because the two bands have different geometry, 



i.e., 3D and Q2D Fermi surfaces. 




FIG. 4: Spatial structure of vortex lattice state for H || ab. 
B — 0.1 and a x /(2a y /v3) = 1-5. The vortex centers are 
located in the middle and at the four corners of the figure. The 
pair potential amplitude [Az,(r)| (a) and |As(r)| (c). Internal 
field 6(r) (e). Zero-energy LDOS N total {E = 0, r) (f) and each 
band contribution N L (E = 0, r) (b) and N S (E = 0, r) (d). 
The peaks are truncated at N = 2iVo. 



IV. VORTEX CORE STRUCTURE FOR 
PARALLEL FIELDS 

We study the vortex core structure and the LDOS 
around the vortex for H \\ ab in order to see how the 
anisotropy of the vortex core shape is affected by 3D or 
Q2D Fermi surfaces. Figure 01 shows the vortex struc- 
ture for H || ab at B = 0.1. In this low field case, 
we use the vortex lattice ratio a x /(2a y /\^3) — 1.5 for 
the unit cell of the vortex lattice. Looking at the am- 
plitude of the selfconsistently calculated pair potential 
|Ai(r)| [(a)] and |As(r)| [(c)] under the given vortex 
lattice ratio, one can notice that the vortex core shape 
is highly anisotropic due to the effect of the Q2D Fermi 
surface. The anisotropy of internal field 6(r) [(e)] is not 
so large. The LDOS on the L-band, N L (E = 0, r) [(b)], is 
highly anisotropic, reflecting anisotropy of | A J c(r)|. How- 
ever, the LDOS on S-band, N S (E = 0, r) [(d)], shows 
isotropic distribution reflecting 3D Fermi surface. The 
LDOS around the vortex core is broad on the S-band, as 
in the case H || c^iLii Since the S-band contribution 
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is dominant at this low field, as discussed in Fig. E{a), 
the total LDOS N tota i(E = 0, r) [(f)] shows the isotropic 
vortex core shape, reflecting N$(E = 0,r). This corre- 
sponds to an almost isotropic vortex core shape observed 
by STMi^4 The highly anisotropic L-band contribution is 
masked by the S-band contribution. Therefore, while the 
pair potential around the vortex core is highly anisotropic 
reflecting the Q2D Fermi surface, the dominant contri- 
bution of the S-band with the 3D Fermi surface gives an 
isotropic vortex core image in the quasiparticle excita- 
tions observed by STM. 




FIG. 5: The same as Fig. 01 but B = 1 and a x /(2a v /V3) = 6. 
It is noted that scales of x and z are anisotropic in this figure 
of a highly distorted hexagonal vortex lattice case at high 
fields. 

Figure [5] shows the vortex structure at a higher field 
B = 1, where we use the vortex lattice ratio of a highly 
distorted hexagonal, a x /{2a y /\/i) = 6. The difference 
of overall views in Figs. 0] and also comes from the 
deformation of the vortex lattice. While the vortex core 
structures of | A L (r) | [(a)], |A s (r)| [(c)] a,ndN L {E = o, r) 
[(b)] are similar to those in Fig. they are seen as if they 
are isotropic in this figure scaled by the coherence length 
ratio 7jj ~ 6. That is, when we see the pair potential 
amplitude and Nl(E — 0, r), the vortex core structure 
is highly anisotropic, reflecting the Q2D Fermi surface, 
both at low and high fields. The anisotropy around the 
vortex in b(r) increases with raising the field, and satu- 
rates to that of 7# as is seen in Fig. Efe). By the effect of 
increasing the field, Ng(E — 0, r) [(d)] is almost flat and 
slightly enhanced at the vortex center. Therefore, the 



vortex core shape in the total LDOS N to tai(E = 0, r) [(f)] 
reflects the spatial structure of N^{E = 0, r) at higher 
fields. 



V. LOCAL SPECTRUM FAR FROM VORTEX 
CORE 

As mentioned above, while both pair potentials |Al| 
and | As | survive up to H C 2, the electronic state in the S- 
band becomes a normal-like state at high fields. The field 
range of this normal- like state is wider for H \\ ab. There- 
fore, it is interesting to see how the superconducting gap 
edges at |Al| and |Ag| in the quasiparticle excitation 
spectrum are smeared by the quasiparticle excitations in 
the vortex states. 
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FIG. 6: Local spectrum ATtotal(-E, r c ) [solid lines] and the 
S-band contribution Ns(E,r c ) [dashed lines] for H \\ ab at 
the midpoint r c of the surrounding four vortices as shown in 
(d). B = 0.01 (a), 0.1 (b) and 1 (c). a x /(2a y /VE) = 1.5 
(B = 0.01, 0.1) and 6 (B = 1). Arrows in the figure indicate 
\Al\ and |As| at each field, presented in Fig. ffi b). 

The local spectrum N tota i(E, r c ) is shown in Fig. HO at 
the point r c far from vortex cores, where the supercon- 
ducting gap structure in the spectrum may survive up to 
higher fields, compared with other points. At a low field 
B = 0.01 [Fig. HJa)], S-band and L-band contributions, 
respectively, have clear gap edges at E ~ As and A^. 
With increasing B [Fig. HJb)], low energy DOS extend- 
ing from the vortex core grows up within the gap edge, 
and the gap edge is smeared. The smearing is eminent 
for the lower excitation gap of As, whose excitation gap 
structure comes from the S-band contribution N$(E, r c ). 
At high field B = 1 [Fig. RJc)], N s {E,r c ) is almost flat, 
i.e., normal-like state. While there are some attempts 
to observe the field dependence of As by STM^^S the 
excitation gap of As becomes difficult to be observed at 
the field range of a normal-like state S-band state. 
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VI. SUMMARY AND DISCUSSIONS 

Based on the quasiclassical theory, the electronic struc- 
ture in the vortex state was studied in a two-band super- 
conductor with the 3D and Q2D Fermi surfaces, mainly 
for H || ab. To demonstrate the appearance of the effect 
due to the geometry of the 3D and Q2D Fermi surfaces, 
we compared the field dependence of zero-energy DOS 
for H || ab and H || c, and analyzed the anisotropy of 
the vortex core structure for H || ab. The 3D and Q2D 
Fermi surfaces, respectively, introduce the isotropic and 
anisotropic contribution, when c- and ab- orientations are 
compared. At low (high) fields, dominant contribution 
to the change of the total DOS is coming from the small 
(large) gap band which has the isotropic 3D (anisotropic 
Q2D) Fermi surface. Therefore, the field-dependence of 
zero-energy DOS shows isotropic (anisotropic) behavior, 
when the field-orientation is changed. This reproduces 
the results of specific heat and thermal conductivity mea- 
surements ,±2i2£i As for the vortex core anisotropy, the pair 
potential around the vortex core is highly anisotropic re- 
flecting the Q2D Fermi surface. However, the dominant 
contribution of the low-energy qausiparticle state, com- 
ing from the S-band with a 3D Fermi surface, gives an 
isotropic vortex core image at low fields. This may be a 



explanation of the almost circular vortex core image by 
STM in MgB 2 even when H |j abM 

As is seen in DOS, LDOS and spectrum, the electronic 
state on the S-band is normal-like state in a wide field 
range at higher fields. However, this does not mean that 
the pair potential |As| on the S-band vanishes in this field 
region. Quite contrary, |As| persists up to H C 2 since su- 
perconductivity in the S-band is induced by he Cooper 
pair transfer from the L-band. The origin of a normal- 
like electronic state on the S-band is the contribution of 
low energy excitations by a supercurrent around vortex 
cores, coupling with the vortex core state of the quasipar- 
ticles. With an increasing magnetic field, the low energy 
excitations are enhanced, and those low energy quasi- 
particles around the vortex are delocalized giving rise to 
large LDOS between vortices. For the small gap in the S- 
band, the low field gives enough excitations to smear the 
gap structure of As, resulting in a normal- like electronic 
state. We demonstrated that MgB2 shows characteristic 
behavior in the anisotropy of c- and afr-axis directions, 
by coupling of the effect of these normal-like electronic 
states in the S-band at high fields and the effect of the 
Fermi surface geometry (3D and Q2D Fermi surfaces) in 
a two-band superconductor. 
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